Actin Isovariant ACT7 Modulates Root Thermomorphogenesis by Altering Intracellular Auxin Homeostasis.

High temperature stress is one of the most threatening abiotic stresses for plants limiting the crop productivity world-wide. Altered developmental responses of plants to moderate-high temperature has been shown to be linked to the intracellular auxin homeostasis regulated by both auxin biosynthesis and transport. Trafficking of the auxin carrier proteins plays a major role in maintaining the cellular auxin homeostasis. The intracellular trafficking largely relies on the cytoskeletal component, actin, which provides track for vesicle movement. Different classes of actin and the isovariants function in regulating various stages of plant development. Although high temperature alters the intracellular trafficking, the role of actin in this process remains obscure. Using isovariant specific vegetative class actin mutants, here we demonstrate that ACTIN 7 (ACT7) isovariant plays an important role in regulating the moderate-high temperature response in Arabidopsis root. Loss of ACT7, but not ACT8 resulted in increased inhibition of root elongation under prolonged moderate-high temperature. Consistently, kinematic analysis revealed a drastic reduction in cell production rate and cell elongation in act7-4 mutant under high temperature. Quantification of actin dynamicity reveals that prolonged moderate-high temperature modulates bundling along with orientation and parallelness of filamentous actin in act7-4 mutant. The hypersensitive response of act7-4 mutant was found to be linked to the altered intracellular auxin distribution, resulted from the reduced abundance of PIN-FORMED PIN1 and PIN2 efflux carriers. Collectively, these results suggest that vegetative class actin isovariant, ACT7 modulates the long-term moderate-high temperature response in Arabidopsis root.

Programmed and self-organized flow of information during morphogenesis.

How the shape of embryos and organs emerges during development is a fundamental question that has fascinated scientists for centuries. Tissue dynamics arise from a small set of cell behaviours, including shape changes, cell contact remodelling, cell migration, cell division and cell extrusion. These behaviours require control over cell mechanics, namely active stresses associated with protrusive, contractile and adhesive forces, and hydrostatic pressure, as well as material properties of cells that dictate how cells respond to active stresses. In this Review, we address how cell mechanics and the associated cell behaviours are robustly organized in space and time during tissue morphogenesis. We first outline how not only gene expression and the resulting biochemical cues, but also mechanics and geometry act as sources of morphogenetic information to ultimately define the time and length scales of the cell behaviours driving morphogenesis. Next, we present two idealized modes of how this information flows - how it is read out and translated into a biological effect - during morphogenesis. The first, akin to a programme, follows deterministic rules and is hierarchical. The second follows the principles of self-organization, which rests on statistical rules characterizing the system's composition and configuration, local interactions and feedback. We discuss the contribution of these two modes to the mechanisms of four very general classes of tissue deformation, namely tissue folding and invagination, tissue flow and extension, tissue hollowing and, finally, tissue branching. Overall, we suggest a conceptual framework for understanding morphogenetic information that encapsulates genetics and biochemistry as well as mechanics and geometry as information modules, and the interplay of deterministic and self-organized mechanisms of their deployment, thereby diverging considerably from the traditional notion that shape is fully encoded and determined by genes.

Amiodarone inhibits arrhythmias in hypertensive rats by improving myocardial biomechanical properties.

The prevalence of arrhythmia in patients with hypertension has gradually attracted widespread attention. However, the relationship between hypertension and arrhythmia still lacks more attention. Herein, we explore the biomechanical mechanism of arrhythmia in hypertensive rats and the effect of amiodarone on biomechanical properties. We applied micro-mechanics and amiodarone to stimulate single ventricular myocytes to compare changes of mechanical parameters and the mechanism was investigated in biomechanics. Then we verified the expression changes of genes and long non-coding RNAs (lncRNAs) related to myocardial mechanics to explore the effect of amiodarone on biomechanical properties. The results found that the stiffness of ventricular myocytes and calcium ion levels in hypertensive rats were significantly increased and amiodarone could alleviate the intracellular calcium response and biomechanical stimulation. In addition, experiments showed spontaneously hypertensive rats were more likely to induce arrhythmia and preoperative amiodarone intervention significantly reduced the occurrence of arrhythmias. Meanwhile, high-throughput sequencing showed the genes and lncRNAs related to myocardial mechanics changed significantly in the spontaneously hypertensive rats that amiodarone was injected. These results strengthen the evidence that hypertension rats are prone to arrhythmia with abnormal myocardial biomechanical properties. Amiodarone effectively inhibit arrhythmia by improving the myocardial biomechanical properties and weakening the sensitivity of mechanical stretch stimulation.

Increased lateral tension is sufficient for epithelial folding in Drosophila.

The folding of epithelial sheets is important for tissues, organs and embryos to attain their proper shapes. Epithelial folding requires subcellular modulations of mechanical forces in cells. Fold formation has mainly been attributed to mechanical force generation at apical cell sides, but several studies indicate a role of mechanical tension at lateral cell sides in this process. However, whether lateral tension increase is sufficient to drive epithelial folding remains unclear. Here, we have used optogenetics to locally increase mechanical force generation at apical, lateral or basal sides of epithelial Drosophila wing disc cells, an important model for studying morphogenesis. We show that optogenetic recruitment of RhoGEF2 to apical, lateral or basal cell sides leads to local accumulation of F-actin and increase in mechanical tension. Increased lateral tension, but not increased apical or basal tension, results in sizeable fold formation. Our results stress the diversification of folding mechanisms between different tissues and highlight the importance of lateral tension increase for epithelial folding.

Stress fiber anisotropy contributes to force-mode dependent chromatin stretching and gene upregulation in living cells.

Living cells and tissues experience various complex modes of forces that are important in physiology and disease. However, how different force modes impact gene expression is elusive. Here we apply local forces of different modes via a magnetic bead bound to the integrins on a cell and quantified cell stiffness, chromatin deformation, and DHFR (dihydrofolate reductase) gene transcription. In-plane stresses result in lower cell stiffness than out-of-plane stresses that lead to bead rolling along the cell long axis (i.e., alignment of actin stress fibers) or at different angles (90° or 45°). However, chromatin stretching and ensuing DHFR gene upregulation by the in-plane mode are similar to those induced by the 45° stress mode. Disrupting stress fibers abolishes differences in cell stiffness, chromatin stretching, and DHFR gene upregulation under different force modes and inhibiting myosin II decreases cell stiffness, chromatin deformation, and gene upregulation. Theoretical modeling using discrete anisotropic stress fibers recapitulates experimental results and reveals underlying mechanisms of force-mode dependence. Our findings suggest that forces impact biological responses of living cells such as gene transcription via previously underappreciated means.

Gastrulation in Drosophila melanogaster: Genetic control, cellular basis and biomechanics.

Gastrulation is generally understood as the morphogenetic processes that result in the spatial organization of the blastomere into the three germ layers, ectoderm, mesoderm and endoderm. This review summarizes our current knowledge of the morphogenetic mechanisms in Drosophila gastrulation. In addition to the events that drive mesoderm invagination and germband elongation, we pay particular attention to other, less well-known mechanisms including midgut invagination, cephalic furrow formation, dorsal fold formation, and mesoderm layer formation. This review covers topics ranging from the identification and functional characterization of developmental and morphogenetic control genes to the analysis of the physical properties of cells and tissues and the control of cell and tissue mechanics of the morphogenetic movements in the gastrula.

Conformational distributions of isolated myosin motor domains encode their mechanochemical properties.

Myosin motor domains perform an extraordinary diversity of biological functions despite sharing a common mechanochemical cycle. Motors are adapted to their function, in part, by tuning the thermodynamics and kinetics of steps in this cycle. However, it remains unclear how sequence encodes these differences, since biochemically distinct motors often have nearly indistinguishable crystal structures. We hypothesized that sequences produce distinct biochemical phenotypes by modulating the relative probabilities of an ensemble of conformations primed for different functional roles. To test this hypothesis, we modeled the distribution of conformations for 12 myosin motor domains by building Markov state models (MSMs) from an unprecedented two milliseconds of all-atom, explicit-solvent molecular dynamics simulations. Comparing motors reveals shifts in the balance between nucleotide-favorable and nucleotide-unfavorable P-loop conformations that predict experimentally measured duty ratios and ADP release rates better than sequence or individual structures. This result demonstrates the power of an ensemble perspective for interrogating sequence-function relationships.

Spatial-omics: Novel approaches to probe cell heterogeneity and extracellular matrix biology.

Complex intercellular interactions as well as biomolecular and biomechanical cues from the extracellular matrix (ECM) profoundly affect cellular functions. Traditional transcriptomic and proteomic approaches have provided insight into disease progression by identifying discrete cellular subpopulations or microenvironmental signatures characteristic of normal or pathological tissues, however these techniques do not examine how a given cellular state relates to its interactions with neighboring cells or its surrounding ECM with multiparametric characterization (i.e. ECM alignment, mechanical forces, crosslinking, etc.). Emerging spatial-omic techniques can provide high-resolution mapping of expression profiles similar to scRNA-seq and mass spectroscopy directly within tissues. The ability to preserve the spatial context of cells within samples, their cellular geometry, as well as their surrounding ECM gives spatial-omics the opportunity to interrogate previously unexplored signaling modalities, which has the potential to revolutionize ECM research and our understanding of fibrotic diseases. In this review, we present current spatial transcriptomic and proteomic techniques and discuss how they may be applied to investigate cell-ECM interactions.

Biomechanical Forces and Atherosclerosis: From Mechanism to Diagnosis and Treatment.

The article provides an overview of current views on the role of biomechanical forces in the pathogenesis of atherosclerosis. The importance of biomechanical forces in maintaining vascular homeostasis is considered. We provide descriptions of mechanosensing and mechanotransduction. The roles of wall shear stress and circumferential wall stress in the initiation, progression and destabilization of atherosclerotic plaque are described. The data on the possibilities of assessing biomechanical factors in clinical practice and the clinical significance of this approach are presented. The article concludes with a discussion on current therapeutic approaches based on the modulation of biomechanical forces.

A Microarray Study of Articular Cartilage in Relation to Obesity and Severity of Knee Osteoarthritis.

OBJECTIVE: To query the transcript-level changes in the medial and lateral tibial plateau cartilage in tandem with obesity in patients with end-stage osteoarthritis (OA). DESIGN: Cartilage was obtained from 23 patients (20 obese [body mass index > 30 kg/m2], 3 overweight [body mass index < 30 kg/m2]) at the time of total knee replacement. Cartilage integrity was assessed using Outerbridge scale, while radiographic changes were scored on preoperative X-rays using Kellgren-Lawrence (K-L) classification. RNA was probed for differentially expressed transcripts between medial and lateral compartments using Affymetrix Gene 2.0 ST Array and validated via real-time polymerase chain reaction. Gene ontology and pathway analyses were also queried. RESULTS: Scoring of cartilage integrity by the Outerbridge scale indicated that the medial and lateral compartments were similar, while scoring by the K-L classification indicated that the medial compartment was more severely damaged than the lateral compartment. We observed a distinct transcript profile with >50% of transcripts unique between medial and lateral compartments. MMP13 and COL2A1 were more highly expressed in medial versus lateral compartment. Polymerase chain reaction confirmed expression of 4 differentially expressed transcripts. Numerous transcripts, biological processes, and pathways were significantly different between overweight and obese patients with a differential response of obesity on medial and lateral compartments. CONCLUSIONS: Our findings support molecular differences between medial and lateral compartments reflective of the greater severity of OA in the medial compartment. The K-L system better reflected the molecular results than did the Outerbridge. Moreover, the molecular effect of obesity was different between the medial and lateral compartments of the same knee plausibly reflecting the molecular effects of differential biomechanical loading.

Impact of single nucleotide polymorphisms in the VEGFR2 gene on endothelial cell activation under non­uniform shear stress.

Single nucleotide polymorphisms (SNPs) in vascular endothelial growth factor receptor 2 (VEGFR2) are associated with coronary artery disease, hypertension and myocardial infarction. However, their association with atherosclerosis remains to be fully elucidated. The purpose of the present study was to determine whether SNPs are involved in atherogenesis, by analyzing their impact on human umbilical vein endothelial cells (HUVECs) under laminar and non­uniform shear stress in a well­established in vitro model that simulates shear stress­induced proatherogenic processes at vessel bifurcations. All experiments were performed using freshly isolated HUVECs. Three SNPs in the VEGFR2 gene (rs1870377 T>A, rs2071559 A>G and rs2305948 C>T) were genotyped and the expression levels of VEGFR2 were semi­quantitatively determined using western blotting. Subsequently, the HUVECs were seeded in bifurcating flow­through cell culture slides and flow (9.6 ml/min) was applied for 19 h, including tumor necrosis factor­α stimulation during the final 2 h of flow. The protein expression levels of VCAM­1, E­selectin and VEGFR2 and the adhesion of THP­1 cells were analyzed in laminar and non­uniform shear stress regions. Data were analyzed for associations with the respective SNPs. The total expression of VEGFR2 was significantly lower under non­uniform shear stress than under laminar shear stress conditions, independent of the genotype. The expression of VEGFR2 between the different shear stress patterns was not significantly altered by the different SNPs. The expression levels of VCAM­1 and E­selectin were lower in the A/A genotype compared with those in other genotypes in rs1870377 T>A and rs2071559 A>G. In conclusion, the results suggested that SNPs within the VEGFR2 gene have a significant impact on shear stress­related endothelial activation.

Cerebellar Control of Reach Kinematics for Endpoint Precision.

The cerebellum is well appreciated to impart speed, smoothness, and precision to skilled movements such as reaching. How these functions are executed by the final output stage of the cerebellum, the cerebellar nuclei, remains unknown. Here, we identify a causal relationship between cerebellar output and mouse reach kinematics and show how that relationship is leveraged endogenously to enhance reach precision. Activity in the anterior interposed nucleus (IntA) was remarkably well aligned to reach endpoint, scaling with the magnitude of limb deceleration. Closed-loop optogenetic modulation of IntA activity, triggered on reach, supported a causal role for this activity in controlling reach velocity in real time. Relating endogenous neural variability to kinematic variability, we found that IntA endpoint activity is adaptively engaged relative to variations in initial reach velocity, supporting endpoint precision. Taken together, these results provide a framework for understanding the physiology and pathophysiology of the intermediate cerebellum during precise skilled movements.

Chronic administration of minoxidil protects elastic fibers and stimulates their neosynthesis with improvement of the aorta mechanics in mice.

Arterial wall elastic fibers, made of 90% elastin, are arranged into elastic lamellae which are responsible for the resilience and elastic properties of the large arteries (aorta and its proximal branches). Elastin is synthesized only in early life and adolescence mainly by the vascular smooth muscles cells (VSMC) through the cross-linking of its soluble precursor, tropoelastin. In normal aging, the elastic fibers become fragmented and the mechanical load is transferred to collagen fibers, which are 100-1000 times stiffer than elastic fibers. Minoxidil, an ATP-dependent K+ channel opener, has been shown to stimulate elastin expression in vitro, and in vivo in the aorta of male aged mice and young adult hypertensive rats. Here, we have studied the effect of a 3-month chronic oral treatment with minoxidil (120 mg/L in drinking water) on the abdominal aorta structure and function in adult (6-month-old) and aged (24-month-old) male and female mice. Our results show that minoxidil treatment preserves elastic lamellae integrity at both ages, which is accompanied by the formation of newly synthesized elastic fibers in aged mice. This leads to a generally decreased pulse pressure and a significant improvement of the arterial biomechanical properties in female mice, which present an increased distensibility and a decreased rigidity of the aorta. Our studies show that minoxidil treatment reversed some of the major adverse effects of arterial aging in mice and could be an interesting anti-arterial aging agent, also potentially usable for female-targeted therapies.

A microphysiological model of the bronchial airways reveals the interplay of mechanical and biochemical signals in bronchospasm.

In asthma, the contraction of the airway smooth muscle and the subsequent decrease in airflow involve a poorly understood set of mechanical and biochemical events. Organ-level and molecular-scale models of the airway are frequently based on purely mechanical or biochemical considerations and do not account for physiological mechanochemical couplings. Here, we present a microphysiological model of the airway that allows for the quantitative analysis of the interactions between mechanical and biochemical signals triggered by compressive stress on epithelial cells. We show that a mechanical stimulus mimicking a bronchospastic challenge triggers the marked contraction and delayed relaxation of airway smooth muscle, and that this is mediated by the discordant expression of cyclooxygenase genes in epithelial cells and regulated by the mechanosensor and transcriptional co-activator Yes-associated protein. A mathematical model of the intercellular feedback interactions recapitulates aspects of obstructive disease of the airways, which include pathognomonic features of severe difficult-to-treat asthma. The microphysiological model could be used to investigate the mechanisms of asthma pathogenesis and to develop therapeutic strategies that disrupt the positive feedback loop that leads to persistent airway constriction.

Meta-Analysis of Genome-Wide Association Studies Identifies Three Loci Associated With Stiffness Index of the Calcaneus.

The stiffness index (SI) from quantitative ultrasound measurements is a good indicator of BMD and may be used to predict the risk of osteoporotic fracture. We conducted a genomewide association study (GWAS) for SI using 7742 individuals from the Taiwan Biobank, followed by a replication study in a Korean population (n = 2955). Approximately 6.1 million SNPs were subjected to association analysis, and SI-associated variants were identified. We further conducted a meta-analysis of Taiwan Biobank significant SNPs with a Korean population-based cohort. Candidate genes were prioritized according to epigenetic annotations, gene ontology, protein-protein interaction, GWAS catalog, and expression quantitative trait loci analyses. Our results revealed seven significant single-nucleotide polymorphisms (SNPs) within three loci: 7q31.31, 17p13.3, and 11q14.2. Conditional analysis showed that three SNPs, rs2536195 (CPED1/WNT16), rs1231207 (SMG6), and rs4944661 (LOC10050636/TMEM135), were the most important signals within these regions. The associations for the three SNPs were confirmed in a UK Biobank estimated BMD GWAS; these three cytobands were replicated successfully after a meta-analysis with a Korean population cohort as well. However, two SNPs were not replicated. After prioritization, we identified two novel genes, RAB15 and FNTB, as strong candidates for association with SI. Our study identified three SI-associated SNPs and two novel SI-related genes. Overall, these results provide further insight into the genetic architecture of osteoporosis. Further studies in larger East Asian populations are needed. © 2019 American Society for Bone and Mineral Research.

Prioritization of PLEC and GRINA as Osteoarthritis Risk Genes Through the Identification and Characterization of Novel Methylation Quantitative Trait Loci.

OBJECTIVE: To identify methylation quantitative trait loci (mQTLs) correlating with osteoarthritis (OA) risk alleles and to undertake mechanistic characterization as a means of target gene prioritization. METHODS: We used genome-wide genotyping and cartilage DNA methylation array data in a discovery screen of novel OA risk loci. This was followed by methylation, gene expression analysis, and genotyping studies in additional cartilage samples, accompanied by in silico analyses. RESULTS: We identified 4 novel OA mQTLs. The most significant mQTL contained 9 CpG sites where methylation correlated with OA risk genotype, with 5 of the CpG sites having P values <1 × 10-10 . The 9 CpG sites reside in an interval of only 7.7 kb within the PLEC gene and form 2 distinct clusters. We were able to prioritize PLEC and the adjacent gene GRINA as independent targets of the OA risk. We identified PLEC and GRINA expression QTLs operating in cartilage, as well as methylation-expression QTLs operating on the 2 genes. GRINA and PLEC also demonstrated differential expression between OA hip and non-OA hip cartilage. CONCLUSION: PLEC encodes plectin, a cytoskeletal protein that maintains tissue integrity by regulating intracellular signaling in response to mechanical stimuli. GRINA encodes the ionotropic glutamate receptor TMBIM3 (transmembrane BAX inhibitor 1 motif-containing protein family member 3), which regulates cell survival. Based on our results, we hypothesize that in a joint predisposed to OA, expression of these genes alters in order to combat aberrant biomechanics, and that this is epigenetically regulated. However, carriage of the OA risk-conferring allele at this locus hinders this response and contributes to disease development.

Effect of mcl-PHA synthesis in flax on plant mechanical properties and cell wall composition.

The high demand for new biomaterials makes synthesis of polyhydroxyalkanoates (PHA) in plants an interesting and desirable achievement. Production of polymers in plants is an example of application of biotechnology for improving the properties of plants, e.g. industrial properties, but it can also provide knowledge about plant physiology and metabolism. The subject of the present study was an industrially important plant: flax, Linum usitatissimum L., of a fibre cultivar (cv Nike). In the study the gene encoding PHA synthase from Pseudomonas aeruginosa, fused to a peroxisomal targeting signal, was expressed in flax plants with the aim of modifying the mechanical properties of plants. Medium-chain-length (mcl) hydroxy acids in flax plants from tissue cultures were detected by GC-FID and FTIR method. The introduced changes did not affect fatty acid content and composition in generated flax plants. Since mcl-PHA are known as elastomers, the mechanical properties of created plants were examined. Modified plants showed increases in the values of all measured parameters (except strain at break evaluated for one modified line). The largest increase was noted for tensile stiffness, which was 2- to 3-fold higher than in wild-type plants. The values estimated for another parameter, Young's modulus, was almost at the same level in generated flax plants, and they were about 2.7-fold higher when compared to unmodified plants. The created plants also exhibited up to about 2.4-fold higher tensile strength. The observed changes were accompanied by alterations in the expression of selected genes, related to cell wall metabolism in line with the highest expression of phaC1 gene. Biochemical data were confirmed by spectroscopic methods, which also revealed that crystallinity index values of cellulose in modified flax plants were increased in comparison to wild-type flax plants and correlated with biomechanical properties of plants.

Selection on the Colombian paso horse's gaits has produced kinematic differences partly explained by the DMRT3 gene.

The Colombian paso horse, the most important horse breed in Colombia, performs specific and particular gaits (paso fino, trocha, and Colombian trot), which display different footfall patterns and stride frequencies. The breed has been selected for gait and conformation for more than 50 years and we hypothesize that this selection has led to kinematic differences of the gaits that can be explained by different genetic variants. Hence, the aims of the study were: 1. To identify if there are any differences in the kinematic and genetic variants between the Colombian paso horse's gaits. 2. To evaluate if and how much the gait differences were explained by the nonsense mutation in the DMRT3 gene and 3. To evaluate these results for selecting and controlling the horses gait performance. To test our hypotheses, kinematic data, microsatellites and DMRT3 genotypes for 187 Colombian paso horses were analyzed. The results indicated that there are significant kinematic and DMRT3 differences between the Colombian paso horse's gaits, and those parameters can be used partially to select and control the horses gait performance. However, the DMRT3 gene does not play a major role in controlling the trocha and the Colombian trot gaits. Therefore, modifying genes likely influence these gaits. This study may serve as a foundation for implementing a genetic selection program in the Colombian paso horse and future gene discovery studies for locomotion pattern in horses.

Genetically determined pattern of left ventricular function in normal and hypertensive hearts.

We sought to assess the inheritance of left ventricular (LV) function using speckle-tracking echocardiography and the impact of hypertension on modifying the genetically determined pattern of contraction in a population of twins. We recruited 92 Caucasian twin pairs, including 74 hypertensive (HTN) siblings. Beyond standard echocardiographic protocol, a speckle-tracking analysis was performed, including global longitudinal strain (GLS). Systolic function, as assessed by ejection fraction, showed moderate heritability (61%); however, GLS showed higher and dominant heritability (75%). Heterogeneity models revealed that there were no differences between the HTN and non-HTN subjects regarding the heritability of GLS. However, the heritability estimates of diastolic function parameters, including early diastolic strain rate, were low. LV systolic biomechanics is highly heritable. GLS shows dominant heritability, despite the presence of early-stage hypertensive heart disease. Early diastolic parameters are rather determined by environmental factors. These findings suggest the presence of a genetic framework that conserves systolic function despite the expression of diastolic dysfunction and may underlie the phenotypic progression towards heart failure with preserved ejection fraction.

Subject-specific Musculoskeletal Model for Studying Bone Strain During Dynamic Motion.

Bone stress injuries are common in sports and military trainings. Repetitive large ground impact forces during training could be the cause. It is essential to determine the effect of high ground impact forces on lower-body bone deformation to better understand the mechanisms of bone stress injuries. Conventional strain gauge measurement has been used to study in vivo tibia deformation. This method is associated with limitations including invasiveness of the procedure, involvement of few human subjects, and limited strain data from small bone surface areas. The current study intends to introduce a novel approach to study tibia bone strain under high impact loading conditions. A subject-specific musculoskeletal model was created to represent a healthy male (19 years, 80 kg, 1,800 mm). A flexible finite element tibia model was created based on a computed tomography (CT) scan of the subject's right tibia. Laboratory motion capture was performed to obtain kinematics and ground reaction forces of drop-landings from different heights (26, 39, 52 cm). Multibody dynamic computer simulations combined with a modal analysis of the flexible tibia were performed to quantify tibia strain during drop-landings. Calculated tibia strain data were in good agreement with previous in vivo studies. It is evident that this non-invasive approach can be applied to study tibia bone strain during high impact activities for a large cohort, which will lead to a better understanding of injury mechanism of tibia stress fractures.